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Outline 


•  Ab-initio  Modeling  of  the  Following  Reactions: 
-O  +  H2  ^OH  +  H 
-O  +  H2O  -^O  +  H2O*,  OH**  +  OH* 

-O  +  HO  ^O  +  HO*,  a  +  OH* 
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Cross  Sections  vs.  Energy  H2(v} 


O  +  H2(0,0)  ->  h  +  oh 


O  +  ->  H  +  OH 
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Comparisons  with  QM 
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O  +  H2(2,0)  ->  H  +  oh 
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O  +  H2(3,0)  ->  H  +  oh 


Adjusting  QCT  to  Match  QM 


O  +  H2(v,0)  ->  h  +  oh 


TDQM  ^App  v(H2)  =  2 
o  QCT  ^App  v{H^)  =  2. 

□  QCT^App  v(H2)=1.75 
A  QCT^App  v(H2)=1.5 
0  QCT^Appv(H2)=1.0 
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Vibrational  Distributions 


QCT  1 8  kcal/mole 


-e-  o  +  H2(v=oj=o)  ->  oh  +  H 
O  +  H2(v=1  J=0)  ->  oh  +  h 

-B-  O  +  H2(v=2J=0)  ->  oh  +  H 
O  +  H2(v=3J=0)  ->  oh  +  H 


Rotational  Distributions 


18  kcal/mole  H2(v,0)->OH(v)=0 
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18  kcal/mole  H2(v,0)->OH(v)=2 


18  kcal/mole  H2(v,0)->OH(v)=1 


(JS/TI  E)  Mp/sp 


Angular  Distributions 


Collision  Energy  18  (kcal/mol) 


o  +  h2(v=oj=o)  ->  oh  +  H 
o  +  H2(v=i  j=o)  ->  oh  +  H 
O  +  H2(v=2J=0)  ->  oh  +  H 
O  +  H2(v=3J=0)  ->  oh  +  H 
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150 


Surprisal  Fits  of  Vibrational 

Distributions 
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H2(v=1 )  +  O  ->  OH(v)  +  H  H2(v=3)  +  O  ->  OH(v)  +  H 


Energy  Fraction  in  OH(v)  Energy  Fraction  in  OH(v) 


Vibrational  Surprisal  Parameter  vs.  Energy 


O  +  H2(v)  ^  oh  +  h 


Overview  of  the  O  +  H2O  System 


•  Principal  channels 

-  O  +  H2O  ^  O  +  H2O* 

-  O  +  H2O  ^  OH**  +  OH* 

•  Experimental  Studies 

-  Shuttle  Observations  (^8  km/s) 

-  Shock  Tube  (~4.3  km/s) 

-  Physical  Sciences  Inc.  (CVF)  measurements  (~8  km/s) 

-  Rate  constants  up  to  ^^2500  K 

-  Recent  measurements 

•  (H20)0'  +  photon  OH*  +  OH**  (transition  state  information) 

•  O  +  HOD(4Voh)  ^  OH*  +  OD* 

•  Theory 

-  Old  'non-reactive'  surface 

•  Johnson  (Quantum) 

•  Redmon  et  al.  (Classical) 

-  Large  disagreements 

•  Measurements  are  sparse  and  Theory  not  really  helpful 


Cross  Section  (cm") 


Calculations  and  Measurements  of 
the  O  +  HjO  Cross  Sections 


O  +  H2O  Reaction  Approach 


•  Short  Term: 

-  Perform  direct  dynamics  caicuiations  with  VENUS 

•  Quick  picture  of  reaction  chemistry  (but  iess  accurate) 

•  AMI,  PM3  ievei  of  theory  (approximate) 


•  Long  Term: 

-  Generate  new  O  +  H2O  surfaces  and  do  reactive 
scattering  caicuiations 

•  Compute  new  high  ievei  ab-initio  reactive  surface 

•  Fit  surface 

•  Run  ciassicai  trajectory  caicuiations  with  VENUS 


*Karkach  and  Osherov, 
JCP,  110  11918  (1999) 


Rate  Constant  (cm  /molec/sec) 


Reactive  Rate  Constant: 


HOD*  +  O  Direct  Dynamics  Resuits: 
Comparisons  to  Experiment 


O  +  HOD(4Voh)  ^  OH*  +  OD* 

Pfeiffer  ef  a/.  J.  Chem.  Phys.  Vol.  113,  7982  (2000):  Experiment 

1)  Excite  OH  bond  in  HOD  with  4  quanta  (about  20  kcai  mol'^  over  barrier) 

2)  Let  excited  HOD  react  with  very  low  energy  (590  cm-^)  O  atom 


3)  LIE  detection  of  product  OH  and  OD 
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•  Theory  yields  somewhat  hotter  distributions  than  measurements 

•  Relaxation  certainly  a  factor  in  experimental  set  up 


Cross  Section  (cm  /molec) 


Consistency  of  Direct  Dynamics 

Approaches 


O  +  H20(0,0,0)  ^  20H 


Direct  Dynamics  Predictions 
for  OH*  Production 


•  Direct  dynamics  cross 
sections  are  comparable  to 
experimental  measurements 
at  km/s 

•  Need  further  analyses, 
especially  spectral 
comparisons  with 
measurements 


20  40  60  80  100  120  140 

Enegy  (kcal/mol) 


state  Specific  Reactive  Cross  Sections 
O  +  H2O  OH(v,J)  +  OH(vVJ  ) 
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Note  that  OH  products  are  formed  in  'hot' 
and  'cold'  populations 
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Cross  Section  (cm  /molec) 


Direct  Dynamics  Results  for 
O  +  H2O  ^  O  +  H2O* 


•  With  classical  approach,  it  is  not  straightforward  to  determine  product 
vibrational  mode 

-  Assume  most  vibrational  energy  goes  into  H20(v2) 

-  Preliminary  results  are  comparable  to  experiment 


O  +  H20(0,0,0)  ^  O  +  H20(0,V2,0) 


Status  of  Direct  Dynamics 

Calculations 


Investigating  possible  H  +  OOH  product  channel  (high 
energies) 

More  calculations  for  better  statistics 
Simulating  spectra  from  final  OH  distributions 

-  Two  temperature  fit 

-  Compare  to  PSI  CVF  measurements  (OH  overtone) 

-  Other  spectral  data 

Continue  comparisons  with  HOD*  +  O  measurements 


O  +  H2O  Potential  Energy  Surface:  S 

Electronic  Structure  Caicuiations 


•  Basis:  TZV  +  2d, If  on  O  +  2p  on  H;  90  total  basis 
functions 

•  lOeSo  CASSCF  +  MP2  (Is, 2s  frozen,  2  lone  pair,  2 
bonding,  2  alpha,  2  anti-bonding) 

•  No  symmetry;  3  state  average 

•  Calculated  AH  =  15.6  kcal  mol'^,  exp.  =  14.7  kcal  mok^ 


Internal  Coordinates  for  Electronic 
Structure  Calculations 


•  As  reaction  proceeds,  x  and  y  change 

-For  Reagents:  x  is  small,  y  is  large 
-For  Products:  x  is  large,  y  is  small 

•  Performed  ^l.e5  fixed  point  calculations 

•  Surface  fitting  code  nearly  complete, 

-  Unit  testing  in  progress 


H 


y  (t)=  out  of  plane  torsion  angle 
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0  =  “cis"  HI -02 


180  =  “trans”  H1-02 


spectator  =  u,q,t 


active  =  x,y,z 


Coordinates  for  “spectator  geometries” 
(parenthesis  is  notation  in  Wu,  Schatz,  et 
ai.  JCP  2000) 


r2(y) 


Electronic  Structure  Calculations:  All 
Coordinates  Fixed  Except  x  and  y 


r1  =  0.957,  a=104.5,  y=0,  (3=150 


state  1  mp2  kcal/mol 
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Intend  to  Follow  Lowest  Three 
Electronic  States 


statel  mp2  n  =  0.957,  a=104.5,  y=0, 
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Surface  Fitting 
Computation  Procedure 


•  First,  compute  electronic  energies  for  selected  geometries 

•  Then,  generate  a  potential-energy  surface  from  data  using  a 
combination  of  pre-determined  functional  behavior  and 
multidimensional  spline  interpolation 

-  Multi-dimensional  spline  routines  from  Princeton  Piasma  Physics 
Lab 

-  Water  and  OH  potentiais  from  Schwenke  et  ai.  [J.  Chem.  Phys., 
1996] 

-  Functional  fits  to  ieast  sensitive  coordinates  from  Schatz  et  ai.  [J. 
Chem.  Phys.,  2000] 

•  Once  complete,  the  potential  corresponding  to  an  arbitrary  set 
of  Cartesian  coordinates  can  be  computed  within  a  classical 
dynamics  simulation. 


Spectator  and  Active 
Coordinates 


The  twelve  Cartesian  coordinates  of  our  four-atom  system  can 
be  uniquely  described  by  a  set  of  six  “internal”  coordinates. 
These  six  coordinates  are  separated  into  three  “active” 
coordinates  (x,y,z),  which  will  be  fit  using  multidimensional 
splines  and  three  spectator  “coordinates”  (u,  0,  x),  which  will  be 
fit  using  pre-determined  functional  forms. 
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0-H  radial  distance  (internal  coordinate  "u")  [a.u.] 


Smoothing  function  (S) 


From  Reactants  to  Products 


The  potential  energy  surface  (“V”)  will  transition  smoothly  from  the 
initial  reactants  state  (“R”)  to  final  products  state  (“P”)  using 
smoothing  functions  (“S”): 


v  =  (\-s,-s,yv^  +  s,-v,+s,-v. 

Reactants  State  Switching  Function  Product  State  Switching  Function 


Potential  Symmetrization 


The  incoming  oxygen  atom  can  bind  to  either  hydrogen  atom  in  the 
water  molecule.  To  model  two  possible  scenarios  using  a  single  set  of 
Cartesian  coordinates,  a  ‘symmeterization’  procedure  is  added  in 
which  the  position  of  the  two  hydrogen  atoms  are  transposed. 
Transitioning  between  the  two  orientations  is  provided  by  a  smoothing 
function  (S). 


Final 


-O, -H,)+0,)+S,-F((ff,-0, -Nj-0,) 


Interim  Fitting  Results  -  Spline  fit  with 

Interpolation 

Z=80° 2=110° 
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Overview  of  O  +  HCI  System 


Many  important  studies,  but  very  little  information  at  high 
velocities 

-  major  channels  O  +  HCI  ^  O  +  HCI*,  Cl  +  OH* 
Experimental  Studies: 

•  Numerous  experimental  studies  of  the  rate  constant  up  to 
~2500  K 

•  Final  state  resolved  measurement  of  O  +  HCI* 

Theory 

-  Large  body  of  work  over  many  years  up  to  the  present 

•  Detailed,  high  level  potential  energy  surfaces  (PES) 

•  Classical  and  Quantum  studies  (rate  constant) 

•  Very  recent  study  (PES+dynamics): 

»  Ramachandran  (Louisiana  Tech),  Peterson,  Bowman 

•  Note:  H  +  OCi  channei  missing  in  aii  studies 

No  clear  picture  of  magnitude  or  relative  importance  of 
channels 


O  +  HCI  Reaction  Approach 


•  Perform  classical  dynamics  calculations 

-  Use  latest  surfaces  provided  by  Professor  B.  "Ramu" 
Ramachandran  (Louisiana  Tech) 

-  Generate  temperature  dependent  rate  constants  to 
compare  with  available  experimental  data 

-  Generate  velocity  dependent,  state-specific  cross 
sections 

-  Analyze  final  state  distributions 

•  Collaboration  with  Ramachandran/Bowman  groups 

-  Quantum/classical  dynamics  study 


Thermal  Rate  Constant  for  0+HCI  Reaction 
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0+HCI  ->  OH+CI  Thermal  Rate  Constants 


1000/T(1/K) 


•  Good  agreement  between 
present  calculations  (QCT) 
and  experiment/QM  results 
for  T<2000K 

-Possible  quantum 
(tunneling)  effects 
< 1500  K 

•  Disagreement  between 
theory  and  experiment  for 
T>2000K  is  not  currently 
understood 

-  0(^D)  contamination 

-  OCI  channel  accessible 
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Velocity  Dependence  of  O+HCI  Reaction 
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O+HCI(v=0,T  =300K)  ->  OH+CI  O+HCI(v=0,T  =300K)  ->  OH(v')+CI 


Velocity  Dependence  of  OH  Vibrational 

Distributions 
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O+HCI(v=0,T  =300K)  ->  OH(v')+CI 


•  Exponential  falloff  of 
vibrational  distributions 
consistent  with  Polanyi 
studies  of  the  effect  of 
barrier  location  of  vibrational 
energy  in  products 

-  Late  barriers  (in  product 
channel)  result  in  mixed 
energy  release  into 
translation  &  rotation 
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Velocity  Dependence  of 
OH  Vibrational/Rotational  Distributions 
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Velocity  Dependence  of 
OH  Vibrational/Rotational  Distributions 


O+HCI(v=0,T  =300K)  ->  OH(v'J')+CI,  =8  km/s 
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0-i-HCI(v=0,T=300K)  ->  OH(v',j')-i-CI,  E,=10  km/s 
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Estimate  of  OH(v’=1)  Rotational 
Temperatures 


O+HCI(v=0,T  =300K)  ->  OH(v'=1  ,T;)+CI  0H(T  )  <->  OH(<j>) 


Based  on  1®*  moment  of  rotational  distributions  (assuming  a 
Boltzmann  distribution  of  OH(v’=1  ,j)),  OH  rotational 
temperatures  approach  6000K 


Fraction  of  Energy 


Energy  Deposition  in  0+HCI  Reaction 


Fraction  of  Energy  in  Various  Product  Modes  Energy  in  Various  OH+CI  Product  Modes 


•  Reactant  translational  energy  is  preferentially  channeled 
into  product  translational  energy 

•  Rotational  excitation  is  substantial  for  U>  4  km/s 
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OH  Differential  Cross  Section 

at  8  km/s 


0+HCI(v=0,T=300K)  ->  OH(v'=0,1 ,2)+CI 
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Relative  Population 


OH(v'=0,l;j')  Rotational  Distributions 

(Averaged  over  Experimental  Distribution) 
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0+HCI(2,1 )  ->  OH(v'=0  j')+CI  0+HCI(2,1 )  ->  OH(v'=1  ,j')+CI 


O+HCL  Energy  Distribution 


0+HCI  Translational  Energy  Distribution 


•  Experimental  velocity  distribution 
of  H radii  et  al.  is  fit  to  sum  of 
Gaussians 

•  Velocity  distribution  is  first 
convolved  with  thermal  motion 
of  NO2  precursor  and  then  with 
thermal  motion  of  HCI  [van  der 
Zande  et  al.,  1991;  Brouard  et  al., 
1992] 

•  Figure  shows  resultant  0+HCI 
relative  energy  distribution 
generated  by  trajectory  code 


E.  (eV) 


HCI  Vibrational  Excitation 


O+HCI(v=0,T  =300K)  ->  HCI(v')+0 


Rotational  excitation  is 
estimated  to  be  5000K- 
8000K  in  the  6-10km/s 
velocity  range 


